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Abstract 

We study the distribution P(u>) of the random variable u = X\/{x\ + x 2 ), 
where x\ and x 2 are the wealths of two individuals selected at random 
from the same tempered Paretian ensemble characterized by the distribution 
ty(x) ~ <p(x)/x 1+a , where a > is the Pareto index and (j){x) is the cut- 
off function. We consider two forms of 4>(x): a bounded function <p{x) = 1 
for L < x < H, and zero otherwise, and a smooth exponential function 
4>(x) = exp(— L/x — x/H). In both cases ty(x) has moments of arbitrary 
order. We show that, for a > 1, P{oj) always has a unimodal form and is 
peaked at uj = 1/2, so that most probably X\ ~ x<i- For < a < 1 we observe 
a more complicated behavior which depends on the value of 5 = L/H. In par- 
ticular, for 5 < 5 C - a certain threshold value - P(oj) has a three-modal (for a 
bounded <j)(x)) and a bimodal M-shape (for an exponential 4>(x)) form which 
signifies that in such ensembles the wealths X\ and x 2 are disproportionately 
different. 

Keywords: Pareto law, Paretian ensemble, Truncated wealth distribution, 
Fluctuations 



1. Introduction 

At the close of the nineteenth century the Italian economist Vilfredo 
Pareto, while studying the statistics of human income and wealth, discovered 
that the distribution \l/(x) of both has a remarkable power-law form, \l/(x) ~ 
A/x 1+a where a > is a parameter which is now called the Pareto 
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index. Thereafter, Pareto's law was verified for various countries, both for 
the wealth of individuals and for their income, and was also observed in 
empirical data for diverse scientific fields (see, e.g., Refs. 0,3 Hi [fl and 
references therein). It was realized that the Pareto index measured from 
the income distribution is typically larger than the index a deduced from 
the distribution of wealth. This is, of course, consistent with the general 
observation that in market economies wealth is more unequally distributed 
than income 0] • For the wealth distribution, the observed Pareto index a is 
as high as 2.3 — 2.5 for developed countries, but may be as low as 0.81 for a 
developing economy like India jsj. 

As a matter of fact, only very few real- world distributions may follow a 
power law over their entire range. As a distribution of wealth, which is large 
but nonetheless finite, the Pareto distribution (PD) is not an exception. In 
the analysis of data, one uses different forms of truncated PDs, such as, e.g., 
the bounded distribution 

T/ N aL a (l/x 1+a , for L < x < H, 

\&(x) = — < - - > ^ 

1 — o a I , otherwise , 

where L and H are lower and upper cut-offs, respectively, and 5 = Lj H < 1. 
On the other hand, some real-world examples are not that abrupt so that 
one seeks to fit the data using a smoother truncation procedure, e.g., 

Hx) = l^i(-L-±), (2) 



2K a (2V5)x 1+a \ x H 

where K a (x) is the modified Bessel function. In both Eqs. (PQ) and (J2J) a 
power-law emerges as an intermediate behavior, so that ^f(x) has moments 
of arbitrary order. 

Before we proceed, it might be useful to remark that the exponentially- 
truncated distribution in Eq. ([2]) is rather ubiquitous and appears in many 
areas in physics. To name but a few we mention the distribution of the first 
passage times for random motion with a bias [3], the distribution of times 
between action potentials (or the ISI distribution) in the integrate-and-fire 



model of neuron dynamics 10j, the distribution of the stopping distances for 



the sliding motion of a solid block on an inclined heterogeneous plane [111 ], the 



avalanche life-time distribution in the mean-field version of the Bak-Sneppen 



model [12J, probability current distribution in disordered systems [13| and 



2 



the distribution of the number of times a particle diffusing in a sphere hits 



its boundary |14|. Thus our subsequent analysis applies to these systems as 
well. 

In this paper we seek an answer to the following question: Suppose in 
a given tempered Paretian ensemble one selects at random two individuals 
with wealths x± and x 2 , respectively How different are X\ and £2? in quest 
for the answer, we introduce a random variable, 

u,= ^^, (3) 
Xi + x 2 

which defines the contribution of one of two individuals to the total wealth 
x\ + X2 of two of them, and calculate its distribution function P(ou). 

Note that the random variables in Eq. ([3]) have apparently first been 



studied in Ref. [15j within a general context of heavy-tailed, non-truncated 
distributions not having a second, nor even a first, moment. It was realized 
that for a < 1 the distribution P(u) has a characteristic [/-shape form 
with a minimum at ou — 1/2. In our language, this means that for such 
unbounded distributions X\ = x 2 is the least probable event, and the wealth 
distribution is highly disproportionate. Note that recently the distribution of 
more general variables of the form xi/(xi+X2 + . . -+xn) have been studied in 
Ref. [l_6| for different heavy-tailed and tempered parental distributions ^(x). 

In our case, all moments of ^(x) exist. Moreover, since (x™) = (x%) 
for arbitrary n, one may generally expect that P(oj) will be a unimodal 
distribution peaked at u = 1/2. We set out to show instead that this is 
not always the case but that surprisingly, the distribution P{uj) exhibits 
a rather rich, sometimes counterintuitive behavior. In particular, we will 
demonstrate that for tempered distributions with < a < 1 the distribution 
P(oo) undergoes a transition from a bell-shaped form (so that x\ = 22 is 
the most probable event) to an M-shaped (for smoothly truncated \l/(x)) 
and a three-modal (for bounded \l/(x)) forms when 5 = L/H becomes less 
than some threshold value 5 C . This furnishes another striking example that 
random variables with tempered heavy tails can be similar, in important 



respects, to random variables with non-truncated heavy tails [17 



2. The distribution -P(^): A general result 

Let $(A) denote the moment generating function of u: 



$(A) = / [ dx 1 dx 2 ^(x 1 )^(x 2 )exp(-X %1 ). (4) 
Jo Jo V X1+X2J 
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Performing the integration over dxi, we formally change the integration vari- 
able X\ — > oj, so that 



f 1 d 



x J x 2 dx 2 ^ (j-r^ x 2j *( Z 2), (5) 

from which we can read off the desired distribution: 

1 f°° oj 

P(cu) = - xdxV{- x)V{x). (6) 

(1-^) Jo 1-w 

One may readily verify that P(oj) is normalized, i.e., duP(u) = 1, once 
is normalized. Note, as well, that for arbitrary ty(x) one finds from 
Eq. © that (w) = J^udcuP(uj) = 1/2, which can be readily seen from the 
following simple argument: 

1 = = <-^) + (-^—) = 2<-^>. (7) 

Xi + X 2 Xi + X 2 Xi + X 2 Xi + x 2 

We will show in what follows that the average behavior is not representative, 
and under certain conditions does not coincide with typical or most probable 
behavior. 



3. The distribution P(u>) for a bounded Pareto law 

Consider first the bounded PD in Eq. (GQ). Note that for such a distri- 
bution, the choice of a random variable as in Eq. ([3]) will automatically lead 
to a symmetric distribution function that depends on L and H only via the 
ratio 5 = L/H, and has a support not on [0,1] but on a smaller interval 
[u c , 1 — uj c ] where, for L < x < H , uj c = 5/(1 + 5). Substituting Eq. ([T]) into 
Eq. (jUJ) we find 



P(u>) 



a 1 



2(1 -5 a ) 2 io 1+a (l 

[~ 2a (u) - 5 2a m- 2a (u) , for cu c < co < 1 - u c , 

(8) 

for oj < lo c or oj > 1 — u c , 



4 



with 



M(w) 



ml a; 



mm 



max 




(9) 



A straightforward analysis shows that for a > 1, regardless of the value of 
<5, P(w) is unimodal with a maximum at u — 1/2. This signifies that in this 
case most probably xi = X2, and hence, both individuals in such a Paretian 
ensemble contribute to the total wealth proportionally. 

The distribution P(oj) for a > 1 is depicted in Figs. [1] and |2j It has a 
cusp-like maximum at u — 1/2; P(oj) becomes narrower and the height of 
the maximum, P(u = 1/2), increases as the Pareto index a increases. On 
the other hand, P{oS) broadens as 5 — > and P(oj = 1/2) decreases. This 
is, of course, quite a plausible behavior which one may expect on intuitive 
grounds. 

When a < 1 the situation appears to be more interesting and complicated. 
As in the previous case, here P{oj) always has a cusp-like maximum at u> = 
1/2. However, there exists a critical value 5 C of the parameter 5 at which 
two inflection points emerge symmetrically in the regions u < 1/2 and ui > 
1/2. For 5 < 5 C each of these inflection points splits into a minimum and 
a maximum so that P(uj) attains a three-modal, jy-shaped form. When 5 
becomes yet smaller than some other critical 5 CC , one observes that the value 
of P{uj) at these local maxima becomes greater than P(u = 1/2). This 
signifies that in such a Paretian ensemble most probably the two individuals, 
selected at random, have disproportionate wealth since X\ 3> x 2 or %\ <i 2 ' 

Let us discuss this case more precisely, focussing on a = 1/2, for which 
the loci of the extrema can be found explicitly. The distribution P(u) for 
this particular case is depicted in Fig. [3j One finds that the critical value of 
S at which two inflection points emerge is 



For 5 < d~ c , each extremum splits into a minimum and a maximum at locations 




(10) 



to 



1 — cu. 



max,2 



1 + 75- VI - 345 + 5 2 
8(1 + 5) 



mm,l 



1 — CO, 



min,2 



1 + 76 + VI - 345 + S 2 
8(1 + <5) 



(11) 
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Further on, at 5 = 5 CC , which is given by 

S cc = ^259 + 144^ - 12^897 + 518^ ~ 0.0108, (12) 

all three maxima become equally high, i.e., P(cu = 1/2) = P(u maX) i) = 
P{u ma x,2) (see the bold black line in Fig. [3]). For yet smaller values of S, 
the maximum at u — 1/2 becomes smaller than the ones at u = u max ,i and 
u = u ma x,2- Hence, in the latter case it is most probable that x\ and X2 
are disproportionately different and either of them, with equal probability, 
dominates the total wealth. 



4. The distribution P{w) for an exponentially-tempered Pareto law 

Consider next the exponentially-truncated ^(x) in Eq. fl2]). Substituting 
W(x) in Eq. (PJD into Eq. ©, we get 

„. , (LH) a 1 f°° dx ( L x 

P[(jJ ) = = — ; rz / — — -exp 



AKl(2y/6) w 1+Q (l - lo) 1 - 01 J x 1+2a \ lox (l-u)H 
! K 2a (2y/5/u(l-u 



(13) 



Note that P(ui) vanishes at the edges u = and u — 1 exponentially fast 
and is symmetric around u — 1/2. Furthermore, one may readily check that 
the first derivative of P(uj) at ou — 1/2 vanishes so that P{uj) is smooth at 
u = 1/2, as compared to the cusp-like behavior observed in case of a bounded 
PD. The question now is whether u = 1/2 is always a maximum. 

Expanding P(u)) in Eq. ffl3|) in a Taylor series around uj = 1/2, we get, 
omitting insignificant numerical factors, 

1 \ 2 / 1 



P(«)~l + » +O^co--rj, (14) 

where the constant before the quadratic term reads 

l-a-2V(5 . 15 

Observe that g is always negative for any a > 1. This means that -P(w) 
is always a bell-shaped function with a maximum at a; = 1/2. Hence, as 
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in the case of a bounded PD, in this case it is most likely that the wealths 
X\ and £2 of two randomly selected individuals will be the same. Typical 
profiles of P(oj) for a = 2 are depicted in Fig. [Hand show that the form of 
the distribution is not very sensitive to the value of S. 

Note next that the case a = 1 appears to be somewhat special since g — > 
when 5—t-O. This means that for sufficiently small values of S the distribution 
P{oj) will become nearly uniform (apart of an exponential truncation in the 
vicinity of the edges) . This trend is clearly seen in Fig. [5] and signifies that 
for such Paretian ensembles the contribution of each of the individuals to the 
total wealth can have any value (inside the edges of the distribution) with 
equal probability. 

We finally turn to the case < a < 1. One can readily see that here, as 
in the case of a bounded PD, there exists a critical value 5 C , defined implicitly 
by the following transcendental equation, 



For 5 > 8 C , P{oj) is always a bell-shaped function with a maximum at oj = 
1/2. For 5 = 8 C , P(oS) ~ 1 except for narrow regions at the edges, where it 
vanishes exponentially. Lastly, for 5 < 5 C , P(oj) has a bimodal M-like shape, 
with maxima close to oj = and oj = 1, oj = 1/2 being the least probable 
value, which is different from the three-modal W^-shaped form observed for a 
bounded PD. This signifies again that for 5 < 5 C and a < 1 the least probable 
event is that x\ = 22 and most likely either of two contributions completely 
dominate the total wealth, i.e., the distribution of wealth is disproportionate. 
In Fig. E] we depict different possible forms of P(oj) for the particular case 
a = 1/2. Note that for an exponentially-truncated distribution the critical 
value of 5 is somewhat larger, 5 C ~ 0.12, than the corresponding value 5 C ~ 
0.0294 found for the bounded PD. 

5. Conclusions 

To conclude, we have studied here the distribution P(oj) of the random 
variable u = X\/(xi+X2), where X\ and Xi are the wealths of two individuals 
selected at random from the same tempered Paretian ensemble characterized 
by truncated Pareto distributions in Eqs. Q or (j2J). We have shown that, 
for a > 1, P{oS) always has a unimodal form and is peaked at oj = 1/2, 




(16) 
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so that most probably X\ ~ Xi- For < a < 1 (which may be observed 
for developing economies, such as, e.g., India jsj) we have encountered a 
more complicated behavior which depends on the value of 5 = L/H. In 
particular, we have realized that for 5 < 5 C - a certain threshold value - 
P{oS) has a three-modal (for a bounded <p(x)) and a bimodal M-shape (for 
an exponential <j)(x)) form which signifies that in such ensembles the wealths 
X\ and Xi are disproportionately different. Such a behavior appears to be 
quite surprising in view of the fact that the parental distributions ^>(x) have 
moments of arbitrary order. Our findings are summarized in the "phase- 
diagrams" presented in Figs. © and (jHJ). 

We finally remark that a similar shape reversal of the distribution func- 



tion has been observed for melting kinetics of a heteropolymer [18j and for 
the Black-Scholes model of the stock options evolution in mathematical fi- 
nance 19j. Both works dealt with the parental distribution ty(x) of the 
^-variables which, in addition to a power-law intermediate tail, has a log- 
normal truncation for large values of x, rather than an exponential one, and 
an exponential one for small values of x. This, of course, yields quantitatively 
different values of the critical parameters but qualitatively the effect is the 
same. We believe that the effect that we have found here is quite universal 
for the tempered Paretian ensembles, and also quite robust. In particular, 



we expect, as it was shown in Ref. [19| for different cut-off functions, that it 



withstands correlations between x\ and Xi, as long as the correlation length 
does not exceed a certain critical value. 
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p(a>) 




Figure 1: P(u) in Eq. © for 5 = 0.001 and different a. 



P(03) 




Figure 2: P(w) in Eq. (|SJ) for a = 1 and different values of (5. 
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P(Q>) 




Figure 4: P(uj) in Eq. flT3j! for a = 2 and different 5. 
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Figure 5: P(ui) in Eq. (|T3| for a = 1 and different 5. 
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Figure 6: P(w) in Eq. (HU) for a = 1/2 and different <5. 
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Figure 7: A "phase-diagram" for the bounded Pareto law. A vertical line a = 1 separates 
the domains in which P(tu) is always unimodal (a > 1) or may have a different shape 
depending on the value of 5. A horizontal line ln(<5 c ) separates the domain in which P{uj) 
is unimodal and the domain in which P(w) attains a three-modal, PF-like shape with the 
maximum at cj = 1/2 being higher than two local maxima near the edges. Above the 
horizontal line ln(5 cc ) one observes a three-modal distribution P(u>) in which two maxima 
near the edges are higher than the one at uj = 1/2. 
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Figure 8: A "phase-diagram" for an exponentially-tempered Pareto law. A vertical line 
a = 1 separates the regimes in which P(ui) is always unimodal (a > 1) or may have a 
different shape depending on the value of 5. A horizontal line ln(<5 c ) separates the domain 
in which P(u>) is unimodal and the domain in which P(uj) attains a bimodal, M-like shape 
with a minimum at uj = 1/2 and two maxima close to the edges. 
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